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Abstract. A simple statistical model for the effects of dephasing on electron transport in one-dimensional 
quantum systems is introduced, which allows to adjust the degree of phase and momentum randomization 
independently. Hence, the model is able to describe the transport in an intermediate regime between 
classical and quantum transport. The model is based on Biittiker's approach using fictitious reservoirs for 
the dephasing effects. However, in contrast to other models, at the fictitious reservoirs complete phase 
randomization is assumed, which effectively divides the system into smaller coherent subsystems, and an 
ensemble average over randomly distributed dephasing reservoirs is calculated. This approach reduces not 
only the computation time but allows also to gain new insight into system properties. In this way, after 
deriving an efficient formula for the disorder-averaged resistance of a tight-binding chain, it is shown that 
the dephasing-driven transition from localized-exponential to ohmic-linear behavior is not affected by the 
degree of momentum randomizing dephasing. 



1 Introduction 

Electron transport through nanosized systems takes place 
in an intermediate regime between classic and quantum 
transport because the characteristic lengths, system size 
L, phase coherence length l^, and mean free path l m have 
the same order of magnitude. By using these lengths, the 
transport can be classified roughly into the following re- 
gimes. A system with phase and momentum randomiza- 
tion ~ l m < L, for which the resistance is ohmic, i.e. it 
increases linearly with the system size. Coherent transport 
through a disordered system l m < 1$ ~ L, for which the 
resistance increases exponentially with the system size. 
A homogeneous system without momentum randomiza- 
tion < l m ~ L, for which the conductor is ballistic, 
i.e. its resistance is length independent. The interplay of 
these completely different transport regimes makes nano- 
sized systems very promising for novel electronic devices. 
Hence the need for a theoretical description, which covers 
all these regimes, is obvious. 

Basically, the nonequilibrium Green's function (NEGF) 
approach for the transport through quantum systems [1- 
3] allows to include arbitrary interactions by self-energies, 
which cause phase and momentum randomization. By us- 
ing the first order self-consistent Born approximation for 
the self-energies, Datta proposed a model [4], which allows 
to adjust the degree of phase and momentum randomiza- 
tion independently However, its application to nanosized 
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systems involves enormous computational efforts. Hence, 
simple models are required to include the effects of de- 
phasing in nanosized systems. 

Biittikcr proposed to use fictitious reservoirs as a de- 
phasing model, where the phase of the electrons is lost 
by absorption and reinjection [5, 6]. These fictitious reser- 
voirs are used nowadays in several models [7-13]. Note, 
that this apparently phenomenological approach of ficti- 
tious reservoirs can be justified from a microscopic theory 
under appropriate approximations [14, 15]. 

Other models include the dephasing effects by stochas- 
tic absorption through an attenuating factor [16] or by 
random phase factors [17, 18]. However, these models are 
still controversially discussed [19] and most of them are 
defined only in the limit of completely momentum ran- 
domizing or completely momentum conserving dephasing. 

In this paper we present a simple statistical model 
for the effects of dephasing on electron transport in one- 
dimensional quantum systems which allows to adjust the 
degree of phase and momentum randomization indepen- 
dently. The model is based on Biittiker's approach but 
in contrast to other models, complete loss of the phase at 
each dephasing reservoir is assumed, and afterwards an en- 
semble average of the quantity of interest (e.g. resistance 
or conductance) over several spatial configurations of de- 
phasing reservoirs is calculated. The model is discussed 
in detail in Section 2 and applied to linear tight-binding 
chains in Section 3. Concluding remarks can be found in 
Section 4. 
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The model is an extension of our previous work [20, 21], 
where the limit of completely momentum randomizing de- 
phasing has been discussed. In this limit the model has 
been used successfully to explain several conductance mea- 
surements of DNA chains [22]. Moreover, by applying it 
to the ID Anderson model with arbitrary uncorrelated 
disorder, the decoherence induced conductivity has been 
studied and even-order generalized Lyapunov exponents 
have been calculated, which were used to approximate the 
localization length [23]. 

2 Model 

We consider a one-dimensional quantum system with the 
single-particle Hamiltonian H. The system is coupled to 
two reservoirs, which act as the source and the drain of 
electrons, see Fig. 1(a). The coherent transport from the 
source to the drain is characterized by a transmission func- 
tion, which is calculated within the NEGF approach [3] 

T SD = 4Tr(lm(S s )Glm{S D )G+), (1) 

where 

G = (E - H - E s - ZdT 1 (2) 

is the Green function of the system. The self-energies Ss/d 
represent the influence of the source respectively the drain 



reservoir on the system. The current is calculated by the 
Landauer formula [3] 

I=~JdET SD (f s -fD) (3) 

where fg/D = f{E — Ms/d) are the Fermi distributions of 
the reservoirs with the electrochemical potential fis/D- 



2.1 Momentum randomizing dephasing 

At first, we briefly summarize our statistical model for the 
effects of momentum randomizing dephasing [20, 21]. 

Dephasing is introduced by a spatial configuration of 
fictitious dephasing reservoirs, where the phase of the elec- 
trons as well as their momentum are completely random- 
ized, see Fig. 1(b). Thus, these dephasing reservoirs divide 
the one-dimensional system into smaller coherent subsys- 
tems with the transmission Tk,k+i- An energy distribution 
function fk is assigned to each of the k — 1, 2 ... Nd de- 
phasing reservoirs, which can be calculated by using the 
conservation of the energy resolved current at the fictitious 
reservoirs 

Tk-l,k {fk — fk-l) — TkM+l (fk+i — fk) = 0. (4) 
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Fig. 1: (Color online) (a) Coherent transport through 
a one-dimensional quantum system. Fermi distributions 
fs/D are assumed for the source and the drain reservoir, 
(b) Dephasing is introduced by fictitious dephasing reser- 
voirs, where the phase and the momentum of the electrons 
are completely randomized. The dephasing reservoirs di- 
vide the system into smaller coherent subsystems and are 
characterized by an energy distribution function fk. (c) 
Momentum conserving dephasing is obtained if two energy 
distribution functions per dephasing reservoir are defined, 
which allow to distinguish between the electrons with pos- 
itive ft and negative momentum ft- 



The solution of this system of Nd linear equations for the 
energy distribution function of the last dephasing reservoir 
is [20, 21] 



fN d 



Nd-l 



=0 



Tk,k+ 



fD 



<L,k= 



(5) 



7fc,fc + i 



This allows the calculation of the current through the sys- 
tem 



I = R - J dET Nd , D (f Nd -f D ) = ^JdE 



fs-h 



Hence, at an infinitesimal bias voltage /is — > the total 
resistance measured in multiples of h/e 2 



N d 



1 



k=0 



T, 



k,k+l 



(7) 



is given by the sum of the subsystem resistances, which are 
determined by the coherent quantum transport between 
neighboring reservoirs. Moreover, we do not only consider 
a single fixed spatial configuration of dephasing reservoirs 
but calculate the ensemble average of the total resistance 
{R} over the dephasing configurations. Note also that by 
assuming an infinitesimal bias voltage the above used con- 
servation of the energy resolved current is equivalent to the 
conservation of the total current. 
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2.2 Momentum conserving dephasing 

Momentum conserving dephasing is obtained in our model, 
if two energy distribution functions per dephasing reser- 
voir are defined, see Fig. 1(c). The and f£~ allow to 
attribute to every electron a definite sign of the momen- 
tum. The absence of energy relaxation then ensures mo- 
mentum conservation during the dephasing process. The 
current conservation constraint causes 2N c i linear equa- 
tions 

— > — y — y <— • 

Tk-l,k.fk-l - TkM+lfk - Rk.k+lfk + Rk-\,kfk = 0) 

Tfc,fc+l/fc+l _ 2fe-l,fc/)t~ - Rk-l,kf£~ + Rfc,k+lfk* = 0) 

(8) 

where Ri_j = 1 — Tij. The energy distribution function of 
the last dephasing reservoir reads 

_ T^ fs + 1 nj^fp ~ N dfD 

JN d ~ Nd ! _ i J J 

and the total resistance at an infinitesimal bias voltage 



R 



N d 

£ 

fe=0 



-N d . 



k,k+i 



(10) 



Hence, for momentum conserving dephasing, the sum of 
the subsystem resistances is reduced by a constant contact 
resistance R c — 1 for each of the fictitious reservoirs. 

If momentum randomizing as well as momentum con- 
serving dephasing exist in the same system, (8) can still be 
used for the momentum conserving reservoirs, if one sets 
fj^ = f£~ = fk for the momentum randomizing dephasing 
reservoirs and uses for these reservoirs a modified version 



of (4) r fc _ 1>fc (f k - fcj 



TkM+l (fk+1 - fk) — 0. 



3 Application to linear tight-binding chains 

In the following we consider linear tight-binding chains of 
N sites with a single energy level per site and coupling 
t between nearest neighbors 



N 



N-X 



(11) 



k=l 



k=l 



To simplify the notation we take t = 1 as our energy unit 
and the lattice spacing a = 1 as our length unit. The 
reservoirs are modeled as semi-infinite chains which cause 
the self-energy (see for example [21]) 



s = 



E 



E 
2 



for E < -2, 
for \E\ < 2, 
for E > +2 



(12) 



on the first and last site of the chain. In the following, 
only the energy band \E\ < 2 is considered because for 
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Fig. 2: (Color online) A dephasing configuration is gen- 
erated by replacing randomly chosen bonds of the tight- 
binding chain with fictitious dephasing reservoirs. 
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Fig. 3: Decision tree for generating dephasing configura- 
tions. 



\E\ > 2 the imaginary part of the self-energy vanishes and 
hence T = 0. The coherent transmission through a chain 
of length N is calculated recursively by using the same 
method as in [23] 



1 



\r N + S 2 s N -i - S(r N -i + s N )\ 



T N (E) 
with the polynomials 



4Im(£) 2 



(13) 



ro = 
r-i 



1, 



0, 



si 
so 



= (E- £j) s j-l 

1, 

0. 



s i-2, 



(14) 



A dephasing configuration is generated by replacing 
each bond of the chain by a dephasing reservoir with prob- 
ability p, see Fig. 2. These dephasing reservoirs are mo- 
mentum randomizing with probability p r and hence mo- 
mentum conserving with probability 1 — p r , see Fig. 3. 

The phase coherence length is defined by the average 
length of the coherent subsystems 



N 



l + (N-l)p' 



(15) 



The average distance between momentum randomizing 
dephasing reservoirs reads 



jdph 



N 



1 + (N - l)pp r 



(16) 



If only of the Nd dephasing reservoirs are assumed 
as momentum conserving, the ensemble average of the re- 
sistance over all dephasing configurations is given by 



-(JV-l)p(l-Pr) 



(17) 



because {N%} = (N - l)p(l - p r 
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The additional resistance due to momentum random- 
izing dephasing 

{AR} = {R( Pr )} - {R{ Pr = 0)} = (N — l) PPr (18) 

increases linearly with the chain length and both dephas- 
ing probabilities. 

Regarding a fixed dephasing configuration as a series of 
incoherently coupled tunneling barrieres, its resistance in 
the limit of completely momentum randomizing dephas- 
ing (7) and completely momentum conserving dephasing 
(10) corresponds to the results in [6, 11, 19]. However, by 
ensemble averaging our statistical dephasing model pro- 
vides a simple formula (18) for the additional resistance 
due to momentum randomizing dephasing, which is also 
valid for partial decoherence p < 1 and partial momentum 
randomization < p r < 1. 



3.1 Homogeneous chains 



Homogeneous chains Sk 
perfect transmission T = 
aged resistance 



= are considered, which show 
1 if li? I < 2. The ensemble aver- 



{R} = {N d + l}-(N-l)p(l- Pr ) 
= 1 + (N - l)pp r = N/l m 



(19) 



is given by the ratio of the chain length and the mean free 
path since in homogeneous chains momentum randomiza- 
tion is caused only by the dephasing reservoirs and hence 
lm — • Obviously, momentum conserving dephasing 
p r = does not cause an additional resistance and retains 
a ballistic conductor {i?,} = 1, whereas momentum ran- 
domizing dephasing p r > leads to an ohmic resistance 
{R} oc N. 



3.2 Disordered chains 

In this section the onsite energies are considered as dis- 
tributed independently according to a probability density 
w(e) with mean value and variance a 2 . 

The onsite disorder causes momentum randomization, 
which for large N can be estimated by the exponential de- 
crease of the coherent transmission T oc exp(— N/X). The 
mean free path l m is then given by summing up the spatial 
rates of the two momentum randomization processes 



1 



1 

idph 



(20) 



The energy resolved transmission of a fixed disorder 
configuration in Fig. 4 shows strong oscillations due to in- 
terference, if the transport is coherent. When dephasing 
is introduced, p > 0, these oscillations vanish because the 
interference between the disordered sites is no longer pos- 
sible. If the dephasing is momentum conserving p r = 0, 
only the oscillations are averaged out, whereas the trans- 
mission itself is reduced, when the dephasing is momen- 
tum randomizing p r > 0. This property of the transmis- 
sion function, which is also reported in [4, 11, 18], clearly 




Fig. 4: (Color online) The coherent transmission through a 
chain of 15 sites with a given disorder configuration shows 
oscillations due to the interference between the disordered 
sites. When dephasing p > is introduced, the interfer- 
ence is reduced and hence the oscillations vanish. If the 
dephasing is momentum conserving p r = only the os- 
cillations are averaged out, whereas the transmission it- 
self is reduced if the dephasing is momentum randomizing 
p r > 0. This clearly indicates the additional resistance. 



indicates the additional resistance due to momentum ran- 
domizing dephasing. 

Now, disorder averages are studied. The disorder-aver- 
aged resistance of a coherent chain of j sites is calculated 
recursively in the same ways as in [23] 



i) s /d ei ...d^n«( e 



P 3 + 2P^ (1 - E 2 ) + Pj_ 3 + E 2 R 3 



4- E 2 



(21) 



with the recursion relations 



Pi = 



(22a) 



/ dei . ..dsj r 2 Y\_w(ei) 
J i=l 
(E 2 + a 2 ) P^x - 2EQ j _ 1 + Pj_ 2 , 

/ dei . . . dej TjTj^i Y[ w(si) 
J i=i 

EPj^-Qj-u (22b) 
/ &£i...&e j r j s ] -x 

J i=l 

= EQj-! - Rj-i - 1 (22c) 

and the initial conditions Pq = 1, P_i = Qo = R\ = 0. 
Note that the disorder-averaged resistance depends only 
on the variance of the probability density w, higher mo- 
ments do not enter. 

If dephasing bonds are introduced with probability p 
in a chain of N sites, the coherent chain appears obviously 
with probability (1 — p) N ~ x . A subsystem with length j < 



R« — 
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Fig. 5: (Color online) Length dependent disorder-averaged 
resistance in the band-center E — (a = 0.5). The resis- 
tance transits from localized-exponential to ohmic-linear 
behavior, if the dephasing probability is sufficiently in- 
creased p > p* > 0. Momentum conserving dephasing 
p r = always reduces the resistance as the fictitious reser- 
voirs do not cause extra resistance but reduce the tun- 
neling distance. Thus, the resistance is minimal if every 
bond is replaced by a momentum conserving dephasing 
reservoir. If the dephasing is momentum randomizing, the 
resistance depends crucially on the chain length, the disor- 
der strength, the dephasing probabilities, and the energy. 



N appears with probability p(l — p) J 1 at one of the two 
chain ends because it is separated by only one dephasing 
bond from the remaining sites. However, inside the chain 
it is separated by two dephasing bonds and appears with 
probability p 2 (l — p) J 1 at one of the N — 1 — j possible 
positions. Hence, the average number of subsystems with 
length j in a chain of N sites is given by 



_ / (i- P r i (2p+(7v-i-j> 2 ) 
^-{(i-p)"- 1 

The averaged resistance reads 



for j < N, 
for j = N. 



N 



({R}) 



3 = 1 



(N-l)p(l- Pr ). 



(23) 



(24) 



Fig. 5 shows the resistance in the band-center E = as 
a function of the chain length TV. Because of the disorder- 
induced (a ~ 0.5) quantum localization, the coherent re- 
sistance p = increases exponentially. However, if the sub- 
system length is sufficiently decreased by dephasing p > 
p* > 0, the resistance transits from localized-exponential 
to ohmic-linear behavior. The critical dephasing probabil- 
ity p* , which has been calculated analytically in our recent 
work [23], can be recognized approximately from Fig. 6 
as the p-value, where the resistance per length diverges. 
Moreover, the critical dephasing probability is indepen- 
dent of the degree of momentum randomizing dephasing 
p r because the linear factor (N— l)p(l — p r ) does not alter 
the fact that the resistance is exponentially increasing. 
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Fig. 6: (Color online) Average resistance per length of a 
chain with N = 1000 sites in the band-center E = 0. The 
critical dephasing probability p* indicates the transition 
from localized-exponential to ohmic-linear behavior. 
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Fig. 7: (Color online) Disorder-averaged transmission of 
the 15-site chain (a = 0.5). Oscillations as in Fig. 4 do 
not appear due to the ensemble average over disorder con- 
figurations. The transmission is always increased by mo- 
mentum conserving dephasing p r — 0, whereas it depends 
crucially on the chain length, the disorder strength, the 
dephasing probabilities, and the energy, if the dephasing 
is momentum randomizing p r > 0, cf. Fig. 5. 



Fig. 5 and the energy resolved transmission in Fig. 7 
show that the disorder-averaged resistance is always re- 
duced by momentum conserving dephasing p r — be- 
cause the fictitious dephasing reservoirs reduce the tunnel- 
ing distance without causing additional resistance. Hence, 
the disorder-averaged resistance is minimal if every bond 
is replaced by a momentum conserving dephasing reser- 
voir. However, if the dephasing is momentum randomiz- 
ing p r > 0, the resistance depends crucially on the chain 
length, the disorder strength, the dephasing probabilities, 
and the energy. Compared to the coherent resistance, the 
decoherent resistance may increase as well as decrease. 
Note that for a fixed disorder configuration, constructive 
interference increases the transmission at some energies 
and hence momentum conserving dephasing reduces the 
transmission, see Fig. 4. 
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4 Conclusions 

This paper introduces a simple statistical model for the 
effects of dephasing on electron transport, which allows to 
adjust the degree of phase and momentum randomization 
independently by using the dephasing probabilities p and 
p r as well as the onsite disorder a. These parameters arc 
related to physical quantities, namely the phase coherence 
length (15) and the mean free path (20). 

The resistance of a tight-binding chain (17) indicates 
that momentum randomizing dephasing causes an addi- 
tive resistance (18), which increases linearly with the chain 
length and both dephasing probabilities. 

Studying a fixed disorder configuration, it has been 
shown in Fig. 4 that only the oscillations in the trans- 
mission are averaged out if the dephasing is momentum 
conserving, whereas the transmission itself is reduced if 
the dephasing is momentum randomizing. 

The disorder-averaged resistance (24) is calculated very 
efficiently by considering the average number of coherent 
subsystems (23) for the dephasing average and a recursion 
formula (21) for the disorder average of the coherent sub- 
systems, which depends only on the variance of the disor- 
der distribution. Using this formula, it has been shown in 
Fig. 6 that the dephasing-driven transition from localized- 
exponential to ohmic-linear behavior is not affected by 
the degree of momentum randomizing dephasing because 
a linear factor does not alter the fact that the resistance is 
increasing exponentially. However, the disorder-averaged 
resistance is always minimal if the dephasing is momen- 
tum conserving because the tunneling distance is reduced 
without causing additional resistance. 

Applying our model to homogeneous chains, it has 
been shown that the resistance is in between the ballis- 
tic and the ohmic regime, whereas for disordered chains, 
the resistance is in between the exponential and ohmic 
regime. Hence, this model covers all three regimes, which 
have been defined in the introduction. 
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